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ABSTRACT 

Supernova 1987A remains the most well-observed and well-studied supernova to 
date. Observations produced excellent broad-band photometric and spectroscopic cov- 
erage over a wide wavelength range at all epochs. Here, we focus on the very early 
spectroscopic observations. Only recently have numerical models been of sufficient de- 
tail to accurately explain the observed spectra. In SN 1987A, good agreement has 
been found between observed and synthetic spectra for day one, but by day four, the 
predicted Balmer lines become much weaker than the observed lines. We present the re- 
sults of work based on a radiation-hydrodynamic model by Blinnikov and collaborators. 
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Synthetic non-LTE spectra generated from this model by the general radiation transfer 
code PHOENIX strongly support the theory that significant mixing of 56 Ni into the outer 
envelope is required to maintain strong Balmer lines. Preliminary results suggest a 
lower limit to the average nickel mass of 1.0 x 10~ 5 solar masses is required above 5000 
km s -1 by day four. PHOENIX models thus have the potential to be a sensitive probe for 
nickel mixing in the outer layers of a supernova. 

Subject headings: line: formation — nuclear reactions, nucleosynthesis, abundances - 
radiative transfer — supernovae: (1987A) 

1. Introduction 

Until the explosion of Supernova 1987A in the LMC, computer models of supernovae were 
usually calculated with assumptions that allowed only the prediction of bolometric light curves. 
The matter temperature was usually assumed to be equal to the radiation temperature, and grey 
opacities were used (Arnett et al. 1989; Blinnikov 1999). Often, the models did not account for 
the effects of radiation on the hydrodynamic gas flow, or the effects of sphericity on the spectrum 
(Schmutz et al. 1990). 

Observations of SN 1987A produced extensive broad-band photometry, especially during its 
earliest phase of expansion. This motivated the development of models to predict its spectral and 
photometric evolution. In our analysis here we make use of the optical data obtained at CTIO 
(Phillips et al. 1988) and the extensive IUE data which was obtained, re-reduced and analyzed by 
Pun et al. (1995). 

While SN 1987A led to the growth of observational data, so too did it lead to an improvement 
in theoretical models. Hoflich (1991) modeled the early phase of SN 1987A, using a pure hydrogen 
atmosphere and accounting for non-LTE (Local Thermodynamic Equilibrium). Eastman & Kir- 
shner (1989) modeled the first ten days of the SN 1987A explosion. H I and He I were modeled 
in NLTE, while metal lines were treated as pure scattering LTE lines. Takeda (1991) used a pure 
H/He envelope in NLTE and attempted to understand the lineshapes of the Balmer lines. 

Our work here utilized a hydrodynamic model developed by Blinnikov and collaborators using 
the STELLA software package (Blinnikov et al. 1998; Blinnikov 1999; Blinnikov et al. 2000). This 
model calculated the light curve for up to six months after the explosion, and includes allowance 
for time-dependent, multi-group radiation hydrodynamics, monochromatic scattering effects, and 
the effects of spectral lines on the opacity (Blinnikov et al. 2000). The hydro dynamical models 
were based on the stellar evolution models of Saio, Nomoto, & Kato (1988), Saio, Kato, Sz Nomoto 
(1988), and Nomoto &: Hashimoto (1988). The procedure for synthesizing SN 1987A spectra from 
Blinnikov et al.'s model and the results are described in § 2. Here we concentrate on the fact that, 
by day 4.52, significant discrepancies begin to appear in the synthetic spectra of SN 1987A with 
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respect to the observed spectra. In § 3 we review earlier works which noted this discrepancy. In § 4 
we show that nickel mixing at early times seems to be the only way to make the synthetic spectra 
agree with observations. In a subsequent paper, we will discuss NLTE effects with regard to the 
results of Blinnikov et al. (2000). 

2. Models 
2.1. PHOENIX 

PHOENIX is a general radiation transfer code that computes temperature, opacity, and level 
populations for each of typically 50-100 radial zones in a moving stellar envelope. PHOENIX allows 
the user to calculate level populations for a multitude of different nuclear species in LTE or NLTE 
(Hauschildt h Baron 1999; Short et al. 1999). For this study, the following species were calculated 
in NLTE: H I, He I, He II, O I, Ne I, Na I, Mg II, Si II, S II, Ca II, Fe I, Fe II, and Fe III. 

The Blinnikov et al. hydrodynamic model consists of an ejecta envelope divided into 300 
zones, with compositions defined for each layer (Blinnikov 1999; Blinnikov et al. 2000). The outer 
layers of this model were re-zoned onto the 75-zone grid that is used by a modified version of 
PHOENIX, preserving the compositions in velocity-space. The inner boundary condition was taken 
to be diffusive, which is well justified in these calculations. The gamma-ray energy deposition rate, 
which is the energy in gamma-rays deposited in the matter in the envelope per unit mass per unit 
time, is calculated for each new zone. The new grid was chosen so as to maintain both composition 
and density profiles with adequate resolution. 

2.2. Spectral Synthesis 

Figures 1-3 show the spectra resulting from PHOENIX computations on days 1.36, 2.67, and 3.59, 
respectively, of SN 1987A. The fits between the IUE spectra and the synthetic UV spectra are very 
good, an indicator of the importance of NLTE effects in the supernova envelope. [Blinnikov et al. 
(2000) noted that the predicted UV flux 5-10 days after shock breakout was much greater than the 
actual flux when computed in LTE.] For days 1.36 and 2.67, the agreement between the synthetic 
and observed optical spectra are also reasonably good. However, by day 3.59 a discrepancy in the 
strengths of the Balmer lines begins to appear in the synthetic spectra, most notably Ha, even 
when NLTE effects are included. Figure 4 shows the PHOENIX-generated optical/near- UV spectrum 
for day 4.52 compared with the observations (Phillips et al. 1988; Pun et al. 1995), and here the 
Balmer lines are clearly much weaker in the synthetic spectrum than they are in the observed 
spectrum. 
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3. Weak Lines 

Schmutz et al. (1990), using various hydrodynamic models, noted the hydrogen atoms appeared 
to be sufficiently excited to create strong lines, despite the fact that the photospheric temperature 
had dropped below the ionization threshold (T « 6000K ). They could offer no explanation for this 
discrepancy. The pure H/He non-LTE model of Takeda (1991) also resulted in weak Balmer lines, 
as did a model by Hauschildt & Ensman (1994), which used an earlier version of PHOENIX with H I, 
He I, Mg II, and Ca II in NLTE and metal line blanketing in LTE. These two groups suggested 
that differences between the theoretical and actual density distributions might be responsible. 

Figure 5 shows temperature vs. velocity for Blinnikov et al.'s models at days 1.36 to 4.52. 
The diamonds on the graph show where is approximately equal to 1 in each model (t^j = 
total continuum optical depth at 5000 A). Assuming the line- forming region lies in the vicinity of 
r std = 1> no * e ^at for days 1 and 2, the gas temperature is well above the ionization threshold 
for hydrogen, thus ensuring the presence of strong Balmer lines. (A good rule of thumb is that 
one sees the strongest lines from the species that is just below the dominant ionization stage. For 
example, Si II lines are strongest if Si III is the dominant stage, etc.) By day 3, the temperature 
is beginning to approach the threshold at r g ^ = 1, and at day 4 the temperature is such that 
sufficient ionization can no longer be assumed in the line-forming region. 

Figure 6 shows the gamma-ray energy deposition predicted by Blinnikov et al.'s hydrodynamic 
model. The total energy deposited per second is shown as a function of r g ^ and the integral is 
performed from the outside proceeding inward. Until we reach the innermost layers of the envelope 
model (t s ^j > 100), the deposition is negligible, and therefore hydrogen atoms cannot become 
sufficiently excited through the absorption of gamma rays. 

4. Nickel Mixing 
4.1. Model Results 

In order for the synthetic Balmer lines to match the strength of the observed lines, either higher 
matter temperatures (to thermally populate the Balmer states) or higher radiation temperatures 
(to radiatively populate the Balmer states in NLTE) would be required. Both of these would 
significantly degrade the quality of the fit, since they would both lead to model spectra that would 
be significantly bluer than the observed spectrum. A third way to strengthen the Balmer lines is 
to populate the states non-thermally, via fast electrons created from the deposition of gamma-rays 
produced by the radioactive decay of 56 Ni. In our models this is controlled by the gamma-ray energy 
deposition rate. In order to obtain a larger energy deposition rate, more mixing of radioactive 56 Ni 
into the supernova's outer envelope than exists in Blinnikov et al.'s model (Blinnikov 1999; Blinnikov 
et al. 2000) is required. In PHOENIX this can be accomplished by replacing the self-consistent energy 
deposition rate from the Blinnikov et al. model with an ad hoc function that simulates a gamma-ray 
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population created by excess nickel mixing. To simulate mixing, we assume a uniform mass fraction 
of 56 Ni distributed homogeneously (i.e., X^j = constant), and local deposition of gamma-rays. The 
new energy deposition rate is therefore much larger (and smoother) than the original (Figure 6), 
allowing for more hydrogen ionization and excitation in the line-forming region. Comparisons 
between the synthetic spectrum and the actual spectrum provide a sensitive probe for how much 
56 Ni must be mixed into the supernova envelope. 

Figure 7 shows the results of modeling the day-4.52 spectrum with constant gamma-ray depo- 
sition, resulting from a nickel mass fraction of Xj^ = 1.0 x 10~ 3 in the envelope. The fit between 
the predicted lines and the actual lines is much better, illustrating the validity of this model. A 
mass fraction of Xj^j = 1.0 x 10~ 3 distributed thoughout our model atmosphere with a mass of 
about 10~ 2 Mq corresponds to a requirement of approximately (1.0 — 1.2) x 10 -5 M of nickel to 
create enough gamma rays above the line forming region of the envelope. The mass is low (and 
corresponds to a lower boundary in velocity space of about 5000 km s _1 ), because at these early 
times the density is high and the continuum optical depth is > 100 already at this mass/velocity 
point, so that spectrum formation occurs in the very outer layers of the ejecta. At later times 
we see deeper due to geometrical dilution. Although the model produced by PHOENIX utilized a 
constant nickel mass fraction (so that the number density of nickel follows the rather steep density 
gradient), this by no means precludes the possibility of inhomogeneous mixing. In the case of that 
the nickel is in clumps, our nickel mass is a lower limit (since the clumps could self-absorb gamma 
rays) and our results indicate that at least some of the clumps must be optically thin to the gamma 
rays that are produced. 

In preliminary lower resolution studies, we found that the synthetic Ha line in day 4.52 was 
slightly redshifted with respect to the actual Ha. This was also observed in synthetic spectra 
for days 3.71 and 5.76 created by the radiative transfer code SYN0W (Jeffery & Branch 1990). 
This indicates that when the line strengths are correct, so is their velocity and demonstrates the 
diagnostic power of detailed spectral synthesis. 

4.2. More Evidence for Nickel Mixing 

The mechanism for the mixing of 56 Ni into the envelope is the subject of intense debate, and 
won't be settled until the core-collapse supernova mechanism is better understood. Kifonidis et al. 
(2000b) and Kifonidis et al. (2000a) suggested that nickel mixing into the hydrogen shell would be 
suppressed in Type II supernovae due to strong deceleration at the He/H boundary. Nevertheless, 
evidence exists to support the results of the PHOENIX models. Early detection of X-rays (Dotani 
et al. 1987) and gamma-rays (Matz et al. 1988), the expansion velocities in infrared line widths 
(Erickson et al. 1988; Witteborn et al. 1989) and in 56 Co gamma-ray lines (Barthelmy et al. 1989), 
the detection of H velocities as low as 800 km s _1 at day 221 (Hoflich 1988), and the late detection 
of the He I 10830 line (Fassia & Meikle 1999) all suggest that significant mixing has occurred in 
the ejecta [see Arnett et al. (1989), McCray (1993) and references therein]. 
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Two-dimensional, axially symmetric models calculated by Pryxell et al. (1991) and by Herant 
& Benz (1991, 1992) suggested that 56 Ni could be mixed into the envelope via finger-like structures 
that are produced shortly after the explosion by instabilities in the matter behind the shock front. 
By t = 90 days, these fingers of nickel had broken up and expanded through radioactive heating 
into low-density bubbles of cobalt or iron, the decay products of 56 Ni. 

Inhomogeneous mixing has already been suggested as an explanation for the so-called "Bochum 
event," an emission satellite on the red side of the Ha line in days 20-100. Utrobin et al. (1995) 
suggested that a high-velocity clump of 56 Ni caused an asymmetric over-excitation of hydrogen in 
the outer envelope. Recent studies by Chandra have detected the presence of diffuse high-velocity 
Fe bubbles at the outer edge of the Cassiopeia A supernova remnant. Hughes et al. (2000) noted 
the detection of Fe bubbles at the outer edge of the supernova ejecta, outside the Si-rich layers, and 
suggested that their diffuse nature could be explained by radioactive heating from decaying nickel. 

Models of the composition of SN 1987A in the nebular phase by Kozma & Fransson (1998a), 
based on the results of various explosion models, used a number fraction of iron of 2 x 10" 5 in the 
hydrogen layer, with a filling factor / « 0.17 — 0.70. This corresponds to a mass fraction on the 
order of 10~ 3 . The results of the simulations of Kozma & Fransson (1998b) indicate that the iron 
mass fraction should be within a factor of 2 of the predicted value. Li et al. (1993) suggest that 
the light curve of SN 1987A during the nebular phase was dominated by ~ 0.078 M Q of 56 Co in 
the form of 50 to 100 clumps. Their models contained nickel clumps that were initially opaque to 
gamma rays with / < 0.01, accounting for « 1% of the total envelope mass. Gamma-ray absorption 
heated the clumps, causing them to expand until / = 0.5. The nickel mass fraction and total mass 
values predicted by the PHOENIX runs are well within the limits imposed by these two studies. 

Detection of the He I 10830 line at times after 10 days suggests that helium is also being 
re-ionized by 56 Ni decay. Fassia & Meikle (1999) applied spectral synthesis models to SN 1987A 
and noted that the best fit to the helium line required at least 3% of the total 56 Ni mass (0.07 M ) 
be mixed to above 3000 km s _1 , but that 56 Ni concentrations should be negligible above ~ 4000 
km s _1 . Fassia et al. (1998) carried out a similar analysis of the Type II supernova 1995V, where 
they noted the best fit to the He line required the 56 Ni density to be constant up to 650-1150 
km s _1 , then fall off by a power law of 8-9 (consistent with our mixing assumption). Only 10~ 6 
M Q of 56 Ni was required above the helium photosphere (« 4250 km s _1 at day 69) in the SN 1995V 
models. Burrows & Van Riper (1995), using Monte Carlo calculations of gamma-ray transport in 
clumpy debris models of SN 1987A's envelope, suggest that up to 50% of the total 56 Ni mass in 
SN 1987A should remain below 1000 km s _1 . 

Herant et al. (1994) suggest that significant iron/nickel inversion may be the result of neutrino- 
driven convection within the fraction of a second just prior to the supernova explosion. Convection 
would destroy the spherical symmetry of the elemental distribution, possibly creating the iron/nickel 
clumps predicted and/or detected in later stages. This theory is supported by 2D hydrodynamic 
simulations performed by Kifonidis et al. (2000a), although in their models the nickel was strongly 
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decelerated at the He/H interface. 

5. Further Research 

Future work with the Blinnikov et al. 1987A model will involve using PHOENIX to synthesize 
the entire time series of the supernova up to around day 100. The detailed time series will indicate 
the location of nickel in the entire envelope, since the spectra probe deeper as time goes on due to 
the geometrical thinning of the envelope. PHOENIX spectra for this period may shed light on the 
later-time phenomena discussed in § 4. Once the whole time series of SN 1987A is modeled, the 
synthetic spectra will be applied to the Spectral-fitting Expanding Atmosphere Method (SEAM) 
in order to calculate the distance to SN 1987A, and hence to the LMC (Baron et al. 1993, 1994, 
1995a, b). Our calculations will also be useful to ascertain whether a 1-D representation of SN 1987A 
is sufficient. So far both breakout models (Blinnikov et al. 2000) and our preliminary calculations 
here indicate this to be the case. 
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Fig. 1. — PHOENIX model spectrum for day 1.36, plotted against observed optical and UV spectra. 
Important optical lines include Ha through H<5, and He I A5876. All optical spectra are taken from 
the CTIO archive (Phillips et al. 1988). All UV spectra are from IUE (Pun et al. 1995). 
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Fig. 2. — PHOENIX model spectrum for day 2.67. Note that He I A5876 has virtually disappeared 
while Ca II H&K lines are beginning to emerge. No observed optical spectrum for day 2.67 (Feb 
26) was available (Phillips et al. 1988). 



- 12 - 





Model, day 3.59 

Optical, day 3.71 
UE, LWP 10210-11, day 3.52 
UE, SWP 30401-02, day 3.53 



2000 



4000 
A (Angstroms) 



6000 



8000 



Fig. 3. — PHOENIX model spectrum for day 3.59. The Ca II lines have become stronger, and there 
is now much greater line blanketing in the UV. Note that the predicted Ha line is now noticably 
weaker than the observed line. 
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Fig. 4. — PHOENIX model spectrum for day 4.52. All synthetic Balmer lines are now much weaker 
than the observed lines and are slightly redshifted. The observed lines are bluer since they formed 
at a higher velocity. 
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Fig. 5. — Temperature structure vs. velocity for days 1.36 through 4.52. The diamonds denote the 
approximate locations of = 1 for each model. 
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Fig. 6. — The total gamma-ray energy deposition rate vs. for day 4.52. The quantity plotted is 
the integrated energy deposited in the envelope starting at the outer edge. The unmixed deposition 
rate is obtained from Blinnikov et al.'s model (Blinnikov 1999; Blinnikov et al. 2000) and the mixed 
rate is that obtained with our assumed mixing. 
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Fig. 7. — PHOENIX model spectrum for day 4.52, with gamma-ray deposition calculated assuming 
local deposition due to a constant nickel mass fraction of 1.0 x 1CT 3 everywhere in the envelope. 



